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Abstract 

It is shown how the boson mapping formalism may be applied as a useful many- 
body tool to solve a fermion problem. This is done in the context of generalized 
Ginocchio models for which we introduce S-, D-, and G-pairs of fermions and sub- 
sequently construct the sdg-boson realizations of the generalized Dyson type. The 
constructed SO(12) and Sp(10) fermion models are solved beyond the explicit symme- 
try limits. Phase transitions to rotational structures are obtained, also in situations 
where there is no underlying SU(3) symmetry. 
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1 Introduction 



It has often been recognized and stated that boson mapping is not only relevant to dis- 
cussions about the relationship between phenomenological boson models and shell model 
type fermion models, but that it also constitutes an attractive many-body formalism in its 
own right (see e.g. Refs. Jl], |2| and references therein). The latter pronouncement is against 
the background that one might profit from the use of boson variables in the description of 
many-fermion problems through the potentially simpler algebraic structures associated with 
bosons. This philosophy has its roots in the very simple observation that bosons commute 
to a c-number, while fermion pairs commute to a non-trivial operator. Correspondingly the 
construction of an orthonormal basis is in general then simpler for bosons than for fermion 
pairs. Except for a few cases (see e.g. Refs. [3-9]), this has not been exploited much as a 
possible simpler route to solution of a many-body problem. 

In this paper we demonstrate how the known algebraic structure associated with s-, d-, 
and g-bosons may indeed be exploited to solve non-trivial fermion models. We analyze and 
discuss an extension of the original Ginocchio SO(8) and Sp(6) models which have also 
been re-interpreted and used in the Fermion Dynamical Symmetry Model (FDSM) ||10|| . 
The extension is based on a closed algebraic structure in terms of fermion pairs which now 
also includes a G-pair and leads to the SO(12) and Sp(10) dynamical symmetries. Such 
an extension has also recently been presented by P. Feng et al. |Ll] who discussed analytic 
expressions for energies in some symmetry limits of the Sp(10) case. Here we are able to 
obtain solutions beyond the symmetry limits by using boson mapping methods. (See also 
Ref. [12 1 for further exploitation of the boson method in this context.) 

Before we introduce and analyze these SDG-generalizations of the Ginocchio (FDSM) 
models, we discuss in Sec.|2| the basic principle underlying the construction of algebraic 
fermion models structured around the concept of favored pairs. We focus namely on the 
respective roles of physical and algebraic angular momenta [fT3|| , thus arriving at a division 
into active and inert parts of the angular momentum from a point of view we regard as 
more fundamental than the usual consideration where such a division is simply introduced 
as a convenience. 

In Sec. |] we then construct generalizations of the Ginocchio (FDSM) models for an 
arbitrary highest angular momentum favored pair and subsequently discuss possible sub- 
group chains and representations (Sec. |j) as well as realizations of these models in the 
nuclear shell-model space (Sec. ||). An equivalent description of collective states in terms of 
the single-particle quadrupole operator is discussed in Sec. || where a comparison with the 
harmonic oscillator quadrupole operator is presented. 

In Sec. we construct the Dyson boson mapping for the pair and multipole operators. 
For the case of highest angular momentum J max =4, namely for SDG-models, we consider 
in Sec. || the role of possible spurious states when the boson analysis is carried out in the 
complete and unrestricted Fock space spanned by boson states of the sdg-type. The possible 
appearance and role of spurious states when the full simplicity of a boson basis is exploited, 
is of course well known and analyzed (see e.g. Refs. 0, [l], |Hf), but there are some detailed 
properties for the present analysis worthwhile to be pointed out and discussed. 
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The results of calculations we are able to carry out are presented in Sec. |9|. In fact they 
can be relatively easily obtained once the fermion problem is transformed to an equivalent 
boson description. We show that the SDG algebraic models include phase transitions from 
anharmonic vibrations to rotations. In particular, we obtain well-pronounced rotational 
solutions in the SO(12) model, which are not related to any SU(3) symmetry limit. In Sec. 



10 we give a summary and formulate some conclusions. 



2 Collective algebraic models and angular momentum 

In algebraic collective models one assumes that collective states (of a given collective mode) 
can be constructed by acting with the generators A of a so-called spectrum generating alge- 
bra [15] A on a reference state |ref). The reference state is a fermion state (not necessarily 
the ground state) and A is an algebra of fermion operators. 

Such models have the advantage of being simple enough to incorporate in them an exact 
conservation of basic symmetries of fermion systems. The simplest way to do that is to 
build the generators of A as irreducible tensors of conserved symmetries. In particular, the 
generators should either conserve or change by a given value the particle number of the 
system, and similarly, either conserve or carry with them a given angular momentum. In 
what follows we concentrate our discussion on the conservation of these two symmetries. 

The conservation of the particle number means that some subset of generators of A 
should commute with the particle number operator and constitute a spectrum generating 
algebra for a given fermion system. Some other generators should decrease or increase the 
fermion number either by one, in which case one unifies the dynamics of even and odd 
nuclei, or by two if one wants to consider only even systems. In this paper we consider only 
spectrum generating algebras built of pair-creation, pair-annihilation and single-particle 
operators. In principle one should consider many-body terms in either of these operators 
which would be a natural extension in the spirit of constructing effective operators resulting 
from a reduction to a small phase space. In the algebraic approach, however, this is a very 
difficult task and only a few developments have been made in this direction (see e.g. Ref. 

H). 

Let us attribute to every generator A NJM of A the angular momentum quantum num- 
ber J and its projection M, and let N enumerate generators with the same J and M, if 
necessary. The assumption that the generators Anjm are irreducible rotational tensors can 
be formulated in terms of the following commutation relation 



[Jim, A njm ,) = -p{J + 1)(JM1M'\JM + M')A NJM+M >, (2.1) 
where Jim are spherical components of the angular momentum operator [17]. We will 



call Jim the physical angular momentum to distinguish it from other angular momenta 
which appear later on. The physical angular momentum is a standard operator acting on 
the coordinates and spins of all particles in the fermion system considered, and it obeys 
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standard so (3) commutation relations^], 

[Jim, Jxm>\ = -V2(1M1M'\1M + M')J lM+M , (2.2) 

We denote this algebra by soj(3) to indicate that it is generated by the operators Jim- 

The commutation relation (|2.1|) is easy to implement whenever the generators A^jm 
are given in terms of the space coordinates and spins. Equivalently, one can construct the 
generators Anjm from fermion creation and annihilation operators by coupling them to a 
given angular momentum. In any case, a commutation relation such as fl2.1|) assures the 
correct transformation properties of generators under rotations of the reference frame. It 
also allows one to determine the angular momentum of a collective state by coupling together 
angular momenta of the generators A^jm- 

Let us assume that in the algebra A there exists a J=l operator A\\m=P\m which acts 



on the labels of generators as an angular momentum [13], i.e. 



[Pim, A NJM >] = -s]J{J + 1)(JM1M'\JM + M')A NJM+M , (2.3) 

Applying this relation to Aum we see that Pim also fulfills the so(3) commutation rela- 
tions as in Eq. ( |2.2| ). We will call P\m the algebraic angular momentum and denote the 
corresponding algebra by sop(3). 

Assumption ( [2 .3D is the key element for constructing collective fermion-pair algebras. 
Even if it constitutes a rather stringent condition imposed upon the algebra A, it is a very 
natural one if one aims at constructing a useful algebraic model. It means that we want to 
include in A the operator which would measure the physical angular momentum. In this 
way the model space, which is generated by acting with generators of A on the reference 
state, will always comprise complete multiplets of the physical angular momentum. By the 
same requirement we also have to assume that P\m is a single-particle operator, similarly 
to Jim- 

Of course, the physical angular momentum Jim itself can belong to the spectrum gen- 
erating algebra A. Eqs. ( |2.1| ) and ( |2.3j ) are then just a duplication of the same property. 
This is the case, for example, in the CM(3) model [0 or the symplectic Sp(3,R) model 
[Pl]|j where generators are constructed from particle coordinates. However, in what follows 
we need not distinguish between the two separate cases when Jim and Pim are equal or 
not, because the former can be considered as a special case of the latter - see also below. 
When Jim and Pim are different operators we have to consider a larger algebra, soj(3)+A, 
comprising both the physical spin J 1M and the spectrum generating generators A NJM - It 
is, however, useful to arrange generators of this algebra in another way. 

Subtracting Eqs. fl2.1p and ( |2.3p we see that all generators A^jm commute with the 



angular momentum operator R\m which is the difference of the physical and algebraic 
angular momentum: 

Rim '■— Jim — P\Mi (2-4) 



*The sign of the right-hand side of Eq. (2.2) corresponds to the standard choice |17| for relative phases 
of components of the angular momentum operator Jim- Similarly, the sign in Eq. ( |2.l| ) corresponds to a 
certain assumed phase convention for the Anjm generators. 
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[RiM>,A NJM ] = 0. (2.5) 
Therefore, we obtain a direct sum A', 

A' = so R {3) © A, (2.6) 

of the sor(3) algebra composed of the operators -Rim, and of the spectrum generating algebra 
A. We will call R\m the inactive or inert angular momentum. 
Moreover, applying Eq. (|2.5|) to Piri=Aum we see that 



[Rim',Pim] — 0, (2.7) 

i.e., in the algebra A' there exist two commuting angular momentum operators. This fact 
has important consequences for the structure of the underlying fermion space. Indeed, 
having at our disposal two commuting single-particle operators we can attribute two angular 
momentum quantum numbers to every single-particle state. More specifically, let us suppose 
that algebra A is built in the spherical shell model space of single-fermion states described by 
the creation operators a^ m .. These states are grouped in multiplets (j-shells) with respect 
to the physical angular momentum and therefore a^ m . are irreducible tensors with respect 
to Jim, 

[Jim, at jmj } = -\fj(j + ^(IMjmjljrrij + M)a+ jm . +M , (2.8) 

where n distinguishes multiplets of the same j. The existence of two commuting angular 
momenta allows one now to group single-particle states in another set of multiplets (r- 
p-shells) each described the by fermion creation operators a^, rm rpmp which are irreducible 
tensors with respect to -Rim and P\m simultaneously, 



[RiM,a+, rmrpmp ] = -y/r(r + l)(lMrm r \rm r + M)a+ Irmr+Mpmp , (2.9a) 



[PlM, K'rmrprrJ = ~ ^ P(P + l )( lM P m P \P m V + M ) a n>rrn r p mp +M- ( 2 - 9b ) 

Here n' distinguishes the r-p-shells with the same values of r and p. Since the physical 
angular momentum is by definition the sum of the active and inactive ones, 

Jim = Pim + -Rim, (2-10) 

(cf. Eq. fl2.4Q ), the j-shells and the r-p-shells are related by the angular momentum cou- 
pling, 

«njm = H ( rm rP m p\j m ) a n'rm r pm p , (2.11a) 
m r m p 



Q'n'rm r pm p 



J2(rm r pm p \jm)a+ jm , (2.11b) 



and the quantum number n comprises three quantum numbers n', r and p. 

We have therefore shown that assumption fl2.3|) induces a very specific structure of the 
fermion space, namely, the grouping of single-particle states into r-p-shells. The values of 
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the quantum numbers r and p are as yet unspecified, but one can easily see that by using 
only the r=0 value one recovers the original structure of j-shells because one then has p=j 
and n'—n. In this case all fermion creation and annihilation operators commute with the 
inactive angular momentum, cf. Eq. ( |2.9aj) . This situation is therefore identical to the one 
where the physical angular momentum Jim itself belongs to the algebra A and the inactive 
angular momentum R\m (P-4 ) vanishes. 

Every r-p-shell contains states with physical angular momenta j—\r—p\, . . . ,r+p. The 
structure of nuclear shells can therefore be embedded in the r-p scheme in several ways. For 
example, the A^-th major oscillator shell contains states with j=\ 1 §,-■-, ° +1 and can 
be identified with the single r-p-shell either for p— Ar ° 2 +1 and r=^ or for p=^- and r= JV ° 2 +1 
The same nuclear shell can also be identified with two or more r-p-shells. For example, if 
one wants to take into account the fact that the highest j= 2Af ^ +1 multiplet of the A^-th 
oscillator shell (the intruder state) is well separated from the other ones, one may consider 
it to be an r-scalar shell and identify the remaining multiplets as another r-p-shell. 

The fact that the fermion space is composed of r-p-shells allows one to construct some 
useful collective spectrum generating algebras A. We may start by constructing a large 
algebra A" in the fermion space of the r-p-shells and then select only those generators 
which are scalar with respect to the inactive angular momentum Rim- In this way we fulfill 
by construction assumption ( j2.3|) . In the next Section we present such a derivation for the 
case of a fermion space composed of a single r-p-shell. A generalization to the case of several 
r-p-shells is straightforward and will not be discussed here. 



The construction in the next Section parallels the original derivation by Ginocchio f20 
who follows the ideas introduced in the pseudo-spin model |21| and uses the coupling of the 
half-integer pseudo-spin i with the integer pseudo-orbital angular momentum k to the total 
angular momentum j. His starting point is the fermion space composed of a single k-i-she\\ 
in which he considers either k or i to be the inactive angular momentum. The cases z=| 
for inactive k and k=l for inactive i are then discussed in detail. In what follows we use 
the notation focused on inactive and active angular momentum r and p, respectively, rather 
then on integer and half-integer angular momenta k and i, as the former classification seems 
to be more fundamental. 



3 Generalization of the Ginocchio model 



In this Section we generalize the Ginocchio prescription |2D| to construct fermion-pair al- 
gebras in such a way that they incorporate fermion pairs up to an arbitrary even angular 
momentum J max . We discuss the application to the case of J max =4, i.e., S, D, and G pairs, 
in some more detail (see also Ref. [11|). 



Consider the fermion Fock space built on a single r-p-shell described by the fermion 
creation operators atmrpm i where r and p are inactive and active angular momenta, respec- 
tively. By the vector coupling of the rm r and pm p indices one obtains the inactive and 



6 



active angular momentum operators in the explicit form 



Rim ■= -fpyV(r +l)/3 ( rm ^m;|lM)(pm p pmp|00)a+ mrpmp a rm /. pm / j , (3.1a) 

m r m' r m p m' p 



Pim ■= -rp\/p(p + 1)/3 ( rm rrml\QQ)(pm p pm' p \lM)a^ mrpmp a rmlrPm/p , (3.1b) 

where 

n — ( — lY~ mr ( — '\] p ~ m Pn (19) 

^rmrprrip \ ± J \ ± J u r,-m r p,-m p ) V "^/ 



and the hat over any angular momentum quantum number is understood as r=\ / 2r + 1. The 
summation over m r and m' r in ( |3.1b| ), and over m p and m' p in ( |3.1a| ) couples to total spins 



R=0 and P=0, respectively. Pim is therefore scalar with respect to the angular momentum 
Rim and vice versa, i.e., Eq. ( |2.7| ) is fulfilled. 

The bifermion creation operators aXn r pmp a rm' pm' , together with their hermitian con- 
jugates a rm > pm ' Orrn r pm p and all commutators of the two sets, give the standard bifermion 
so(2(2r + l)(2p + 1)) algebra. The spectrum generating algebra A can be obtained as a sub- 
algebra of so(2(2r + l)(2p + 1)) composed of combinations of generators which are scalars 
with respect to the inactive angular momentum Rim- Therefore, A is given by the bifermion 
creation operators, 

F jm = J\ (rm r rm' r \00)(pm pP m' p \ JM)a+ mrpmp a+ <pm|) , (3.3) 

m r m' r m p m' p 

their hermitian conjugates -Fjm, and all their commutators. In expression (|3.3j ) we have 
included an arbitrary normalization factor g in terms of which the two-particle states are 
normalized as (cf. Ref. [|TJ) 

(0\F JM F+, M ,\0) = g8jj,5 M M>. (3.4) 

Since the physical single-particle (half-integer) angular momentum j is the sum of the 
inactive angular momentum r and of the algebraic angular momentum p, we only consider 
here cases when either r is integer and p half-integer, or vice versa. 

As shown in the Appendix, the algebraic structure of pair operators Fj~ M is equivalent 
to that of fermion pairs occupying a single-p shell for half-integer p, or to that of boson pairs 
occupying a single-p shell for integer p. The former gives the orthogonal algebra so(4j9+2) 
while the latter gives the symplectic algebra sp(4p+2). The single-p shells are repeated 2r+l 
or (2r+l)/2 times, respectively, and the copies are numbered by the magnetic quantum 
number m r of the inactive angular momentum r. 

In a single-p shell for either half-integer or integer p, bosons or pairs of fermions can 
only be coupled to even total spins J. Therefore, the spectrum generating algebra A is in 
both cases comprised of pairs of fermions F]~ M for J=0,2,. . . , J max , where J ma , x =2[p} (doubled 
integer part of p) . 

Algebraic properties of the models constructed here are, therefore, the same as those of 
particles restricted to a single-p shell. On the other hand, the Hamiltonian is a scalar with 
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respect to the physical angular momentum and couples different copies of the single-p shells, 
giving rise to rather rich dynamic properties. 

Commutation relations of A can be obtained by either considering the underlying single- 
p structure or by an explicit calculation. The latter method can be facilitated by using 
useful techniques developed in Ref. p2| . One then obtains the so(4p+2) and sp(4p+2) 
commutation relations for half-integer and integer p, respectively, in the form: 



J2M21 



^ JxMiJ[M[ r J[M^ 

J[M[ 



where the structure constants read: 



C 



J 2 M 2 J' 2 M! 2 
JiM!J[M[ 



%J 1 kJ[J , 2 £ {J\M x J[M' x \J h ' M"){J 2 M 2 J' 2 M' 2 \J h 'M" 



r 



J"=0,l,...,2p 



(3.5) 



V V Ji 
p p J[ 

J 2 J 2 J" 



r 



— (pm p pn p \J 1 M 1 )(pm p pn p \J[M[)(pm' p pm p \J 2 M 2 )(pn' p pn p \J^ 

ip Trip Tip 

(3.6) 

These same commutation relations can be presented in a different form by explicitly 
introducing the one-body operators of A, 



Pjm = h ( rm rrm' r \00)(pm p pm' p \JM)a^ mrpmp a rmlrPm/p , 

m r m' r mpm' p 



(3.7) 



which can have angular momenta J— 0,1,2,. . . ,2p. Again, we include in the definition of 
Pjm an arbitrary normalization factor h, which can facilitate the comparison of different 
conventions used in the literature. For example, the factor which give s the operator P\m 
normalized as an angular momentum, Eq. ( |3.1b| ), is h=—fp^Jp{p + l)/3. The so(4p+2) and 
sp(4p+2) commutation relations now read 



7 jm,Fj, m , 



h f 



x 



]T (J,-MJ'M'\J"M")l 

0,1,.. .,2p { 



J J' J" 
P P P 



P 



J"M" 



(3.8a) 



D JM, Fj'M' 



-2h(-l 



|2p 



JJ' 



X 



{JMJ'M'\J"M") j 

7"=0,2,..,2[p] I 



J J' J" 
P P P 



r J"M" 



(3.8b) 



[Pjm, P. 



J'M' 



h(-l 



l 2 P 



JJ' 



X 



£ ((-1 

J"=0,l,...,2p 



,J+J' 



){JMJ'M'\J"M"^ 



J J' J" 
p p p 



Pj"M" 

(3.8c) 
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We conclude this section by presenting the generators A in terms of physical angular 
momenta of single particle states, i.e., in terms of creation and annihilation operators of the 
j-shells, Eq. ( J2.9ap . Recoupling the angular momenta one obtains the pair operators 
and the one-body operators ( |3.7|) in the form: 



Pjm = ^(-l)^^M ( {i * l)(aiaJ J ' (3.9b) 



r 

JIJ2 



p p r J V - ,1 32 J M 



4 Subgroup chains and representations in the model 
space 



From the commutation relations ( p.8|) one can deduce the subgroup chain structure of the 
generalized Ginocchio (FDSM) model. For the case K = 0, i = p one has 

SO(2(2i + 1)) D XJ(2t + 1)D Sp(2? + 1) D SO(3) (4.1) 

and 

SO(2(2i + 1)) D Sp(2z + 1) <g> SU(2) D SO(3) ® SU(2) , (4.2) 
while for the J = 0, p = k case the subgroup chain structure is 

Sp(2(2A; + 1)) D XJ{2k + 1) D SO{2k + 1) D SO(3) (4.3) 

and 

Sp(2(2A; + 1)) D SO(2k + 1) ® SU(2) D SO(3) <g> SU(2) (4.4) 

The groups U(2p+1) in relations (|4~T| ) and ( |4.3|) are generated by the multipole operators 
Pjm in (|3.7| ). The generators of the groups Sp(2p + 1) and SO(2p + 1) are the multipole 
operators Pjm with J odd and p = z in the former and p = k in the latter case. The angular 
momentum group SO(3) is generated by P\m an d the seniority group SU(2) in the chains 
d3) and (0) by F+, F 0Q and P o- 

For the SO (8) model another algebraic chain is possible, namely one which includes 
SO(7). This possibility is linked to the accidental zero of a 6j-symbol, 

= (4.5) 




and the corresponding simplification of the relations ( 3.8|) . 

We now turn to the SO(12) and Sp(10) models in some more detail. For these models 
the fermion model space is obtained by successive action of the respective S, D, and G 
pair creation operators on the fermion vacuum. These states could also be obtained by 
combining configurations of the active and inert spaces. Only those states of the inert space 
are considered in which single-particle pseudo-orbital angular momenta or pseudospins are 
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pairwise coupled to 0. For the k = 2 and i — | active parts, the basic building blocks are 
thus characterized by the symmetries [11] and [2] in the respective inert spaces. The total 
wave function should of course be antisymmetric and if states are classified according to 
Young tables, conjugate Young schemes in the active and inert spaces are combined. As the 
wave function in the inert space has either I or K = 0, while in the decomposition of the 
product of two conjugate Young schemes the totally antisymmetric representation occurs 
just once, states in the fermion model space are uniquely denoted by the allowed labels 
corresponding to the active space. 

For the Sp(10) case, the k = 2 active part lies in the SU(5) space whereas the inert part 
is in the SU(hP) space. The possible representations for the inert part are given by the 
Young labels [fi, / 2 , . . . , / 2 n/ 5 ] with 5 > f x = f 2 > / 3 = U > • • • > /an/s-i = /mi/5 > and 
fx + /2 + • • • + /2fi/5 = 2iV, iV denoting the number of fermion pairs. Allowed conjugate 
representations in the active SU(5) space have the Young labels [2N — 2p — Aq — 6r — 8s, 2p + 
2q + 2r + 2s, 2q + 2r + 2s, 2r + 2s, 2s] = [2N - 2p - Aq - 6r - 10s, 2p + 2q + 2r, 2q + 2r, 2r\. 
Here, p, q, r, and s are non-negative integers such that Ap + 6q + 8r + 10s < 2N. Moreover, 
the dimensionality of the inert space restricts SU(5) representations in the active space to 
2N — 2p — Aq — 6r — 8s < ^P. This is just the Pauli principle restriction. 

Analogous considerations hold for the SO(12) case with the i = | active part in the 
SU(6) space. The possible representations for the SU(^P) inert part have the Young labels 
[2/ l5 2f 2 , 2/ 2Q/6 ] with 6 > 2ft > 2/ 2 > . . . > 2/ 2Q/6 > and 2f x + 2/ 2 + . . . + 2f 2n/6 = 
2N. Conjugate representations in the active SU(6) space are denoted by the Young labels 
[N - j - 2k,N - j - 2k,j + k,j + k,k,k] = [N -j -3k,N -j -3k,j,j,0}. Here, j and 
k are non-negative integers such that 2j + 3k < N. The dimensionality of the inert space 
now restricts SU(6) representations in the active space to iV — j — 2k < ^P. 



5 Links to the nuclear shell-model space 

The major nuclear shells consist of normal-parity subshells, with values of single-particle 
angular momentum ranging from j=\ to j= 2Af ° -1 where A^ is the radial quantum number, 
and of the unique-parity orbitals with j= 2jV ° +3 . In Table p] we give for iVo=3, 4, 5, and 6 
possible identifications of the normal-parity subshells with the r-p-shells. All cases involving 
fermion pairs up to J max =6 are enumerated. In addition, we also list identifications which 
leave apart the smallest subshell of j=\- This would require treating this subshell as an 
r-scalar. In all cases the unique-parity orbitals have to be treated separately and in the 
algebraic models are always considered to be the r-scalars. Identifications exploited in the 



FDSM |T(J are denoted by asterisks. 

One can see that the Sp(10) SDG fermion-pair algebraic models can be constructed for 
N Q =A and 5, and the SO(12) SDG models for A^ =5 and 6. It is interesting to note that 
both models can be constructed in the same N —5 shell-model space. If one leaves out the 
2=\ shell, the SO(12) model can also be constructed for N =A and the Sp(10) model can 
be constructed for N =6. 

The SD SO (8) and Sp(6) models of the FDSM can therefore be generalized to incorporate 
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the G fermion pairs. In particular, in the Nq=4, 5, and 6 major shells all four models, SO (8), 
Sp(6), SO(12), and Sp(10), can be constructed and compared, provided that in some cases 
one accepts to leave apart the smallest subshell. In Sec. § we discuss the Nq=6 case in some 
more detail. 



6 Quadrupole single-particle operator 

In the algebraic collective models discussed all possible iV-pair states can be obtained by 
acting on the fermion vacuum with the pair creation operators, 

\9)=F+ Ml ...F% Utf \0). (6.1) 

Therefore, one often considers these operators as basic building blocks which allow for a 
complete description of the collective space. On the other hand, as proposed by the Single- 
particle Coherent Excitation Model (SCEM) |2"3J1 , one can also consider the quadrupole 
single-particle operator Qm=Pim and the monopole pair S + =F^ as another set of build- 
ing blocks. Then any other pair creation operator F]~ M can be expressed as a multiple 
commutator of Qm with S + (see Eq. (|3.8|b)), 

F? M ~([Q-[Q-...[Q;S + ) ...]]) jm , (6.2) 

with an appropriate angular momentum coupling. In particular, the stretched coupling 
scheme can be used to raise the angular momentum of fermion pairs: 

(»■ F f€Z = -2V5M-lf>! { I J p ' ; 2 } F} +%M+M ,. (6.3) 
In this way, all pairs for J>0 can be obtained because the 6-j symbol involved here does not 



vanish |T7j as long as J<2[p\. The complete collective algebra can therefore be recovered 
from the operators Q and S + . 

As a direct consequence of this fact we observe that all collective states can be obtained 

as 

I* > = Q Ml ■ ■ ■ Qm Ki S + . . . Q Ml ■ ■ ■ Qm Kn S + \0), (6.4) 

i.e., as a sequence of N monopole-pair creation operators S + interspersed with quadrupole 
operators Qm, and the sets of states ( |6.1|) and (|6.4|) are identical. The total number of 
quadrupole operators in the sequence can be large, and is only limited by the angular 
momentum accessible in the collective space. 

Similarly, all generators Pjm constituting the core subalgebra can be obtained from 
multiple commutators of Qm, see Eq. (|3~8|c). Hence, in cases when the 5* + -pair condensate 
contains components in all irreducible representations of the core subalgebra, all collective 
states can also be obtained from 

\V) = Qm 1 ...Qm k (S + ) N \U)- (6-5) 
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For the SD algebraic models we have at our disposal the set of six building blocks S + and 
D^ I =F2 M and another equivalent set of six building blocks S + and Qm- However, for models 
which use more pairs, like the SDG SO(12) and Sp(10) models discussed in the present paper, 
the number of pair creation operators Fj M , which is equal to (J m ax+l)(</max+2)/2, increases 
very rapidly with the maximum angular momentum J max . On the other hand, these same 
six operators S + and Qm are always sufficient to obtain all states for arbitrary value of J max - 
The set of building blocks S + and Qm allows therefore for a more economic description of 
the collective space. 

In particular, we may compare different algebraic models by comparing properties of 
the corresponding quadrupole operators Q. In fact, in the algebraic models discussed the 
monopole pairs are always equal to the seniority-zero pairs (cf. Eq. fl3.9a|) ), 

S+ = \/f ~ H a im & Jn> ( 6 - 6 ) 

with amplitudes independent of j. Therefore they are uniquely determined if the set of 
j-shells is fixed. On the other hand, in the given space different models have different 
quadrupole operators Qm, 

QM = hJ2 Cji ,h ( a h ~ a h ) M , ( 6 • 7 ) 

3132 



defined in terms of structure constants Cj u j 2 , see Eq. ( p.9b|) . Since j\ and j 2 can differ by 



2 at most, Q is determined be the set of amplitudes Cjj, Cj-ij, and Cj_ 2 j- 

For N =6, and for the orthogonal and symplectic models listed in Table | these ampli- 
tudes are presented in Fig. |]. Each of the six algebraic models which can be defined in 
this shell- model space is completely determined by 14 amplitudes (or 12 if the 3=\ orbital 
is excluded) shown in the Figure. They represent the set of unique values required by the 
algebra closure conditions (|3.8|) . 

In principle, the collective space in this shell-model space should be determined by dy- 
namic conditions starting from realistic two-body interactions and using an appropriate 
effective-operator theory. Such procedure would in general give quadrupole operators more 
or less departing from the algebraic solutions and could either provide these models with a 
microscopic justification or disprove. This programme has not yet been attempted, prob- 
ably because of a relative simplicity of the algebraic models, which allows one to easily 
obtain theoretical predictions for experimental quantities, whereas the difficult question of 
microscopic justification would require a much larger effort. 

In lack of microscopic derivations of the quadrupole operator, we compare the algebraic 
quadrupole operators with the mass quadrupole- moment operator. The amplitudes C Jlj2 
of the v 2 Y2m{0, 4>) operator reduced to one spherical harmonic-oscillator shell read 

C jllh = - J=/ W2 (20i 2 0|Z l0 )(-l)^l^^i | £ ^ I } , (6.8) 
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where Ij 1; j 2 is the radial integral, 

2 " 7 " dri2 ^W r ^ (r) - ^ j -^N -l + l)(N + l + 2) for |/ 1 -/ 2 |=2 ' 

(6.9) 

and l=(li+l2)/2. Since the formula concerns one oscillator shell, the orbital angular mo- 
mentum / has a unique value for every value of j. Note that for li=l 2 the formula gives 
the root mean squared radii (r 2 ) of the harmonic-oscillator states, independent of I. The 
harmonic-oscillator amplitudes (|6.8| ) are plotted in Fig. [I] in units of this value of (r 2 ). 

For the SD SO(8) and Sp(6) models the diagonal amplitudes Cjj differ very much from 
the harmonic oscillator values. However, when the G, and then I pairs are added, the 
agreement becomes better. The same observation also holds for the amplitudes Cj-ij. 
Here one may in principle change the sign of every amplitude by an appropriate changes of 
phases of single-particle states. For the SDGI models, and also for the harmonic oscillator, 
the phases in Eqs. ( |3.9b| ) and (|6.8|) are such that these amplitudes are positive. For the 
other models one could have fixed the phases in a different way and change signs of Cj-ij, 
but even then the agreement with the harmonic oscillator values would not become much 
better. For the algebraic models the signs of amplitudes Cj-2j are different than in the 
harmonic oscillator. 

Of course the harmonic-oscillator amplitudes are not necessarily the best possible dy- 
namical values for the quadrupole moment operator in a given major shell. However, this 
operator has often been used in the conventional shell model and is known to provide a 
fair first approximation for quadrupole collective phenomena. Our results indicate that the 
algebraic models restricted to S and D pairs are rather remote from the harmonic oscilla- 
tor estimates. The desirability to include higher angular momentum pairs has also been 



suggested by many other dynamical analyses, e.g. Refs. 25, 26 



7 Boson realization of the model 

The Dyson boson realization of an arbitrary fermion algebra composed of pair creation 
operators, their hermitian conjugates, and of all independent commutators has been given 
in Ref. [|]]. For the structure constants of the so(4p+2) and sp(4p+2) algebras given in Eq. 
( |3.6|) we obtain the following mapping: 

2a ( P P J ) 

J M J ' { J 2 J x J 1 J 

(7.1a) 

F JM <— > B JM , (7.1b) 

— y(-iy +2p+i i:^{ J ; J ; ^A)^ (7.1c) 

where B JM =(-l) J - M Bj_ M . 



13 



The boson mapping of the so(4p+2) and sp(4p+2) algebras has the same structure as 
the generalized Dyson boson mapping || 27[] of bifermion operators (the notation is the 
same as in Ref. [§]) 



a a aP 



apa a 
a a ap 



B a9 B 



B ap _ B *e B PP B 



6p, 



(7.2a) 
(7.2b) 

m . (7.2c) 

In this case, each of the single-particle indices a, /3, . . . denotes all quantum numbers of the 
fermion creation and annihilation operators, i.e., either rm r pm p , or in the vector coupled 
form, Eq. ( |2.11a| ), rpjm. The mapping ( \l.2\ ) can be transformed to the collective variables 
by introducing collective amplitudes relevant to our model, 



y rpj 



V2 



--RPjih 




which fulfill the orthogonality and completeness relations 

,RPJ^,R'P'J sr A 

" Orr/0 



ji <h 



>RR n) PP' 



E RP.J^,RPJ 



RP 



l+Si 



(7.3) 

(7.4a) 
(7.4b) 



The mapping of the spectrum generating algebra A, which corresponds to R — 0, can 
be achieved simply by dropping all other collective amplitudes with R ^ due to the 
skeletonization theorem pR [29| , and one obtains Eqs. (ITT 



Alternatively, we can construct this boson mapping making use of generalized coherent 
state approach [0O| |. It is based on a one-to-one correspondence between the coherent states 
in the fermion space and in the boson space 

\C) = exp^C^F+JlO) — \C) = e X p(E^M^k)iO)- 



The action of boson operators on the coherent state is simply 

d 



B\ M \C) 



dC 



-|C), 



JM 



B JM \C) = C* JM \C). 



Similarly the fermion collective pair Fj~ M acts as 



Ff M \C) 







9C* JM 



\C). 



(7.5) 

(7.6a) 
(7.6b) 

(7.7) 



On the other hand the action of the other fermion operators on the fermion coherent state 
has to be evaluated using the commutation relations ( p.8|) . Any fermion state and the 
corresponding boson state |\I/) can be represented using the coherent states as (C\^/) and 
(C|\I/), respectively. Then we can find boson images of fermion-pair operators identifying 

C JM <— B\ M , -J!— <— ^ B JM (7.8) 



dC 



JM 



which results in ( |7.1| ). 
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8 Spurious states 



Bosonization of any fermion model from expressions ( Q|) naturally leads to the identification 
of a corresponding boson Fock space. While the so-called physical states, obtained by 
repeated action of the images (|7.1]a) on the boson vacuum, still obey restrictions due to 
the Pauli principle (as discussed in Sec. f|), this does not hold for general states in the 
boson space and in particular not for typical boson basis states of the complete Fock space. 
Diagonalization in this basis therefore generally leads to the occurrence of spurious states. 
This situation and the identification of spurious states is discussed at length in Refs. || [[]. 



14, a 31 



For the SO(12) and Sp(10) models the associated ideal boson space spans the symmetric 
representation [N] of the sdg boson U(15) group. In the U(15)dSU(5) reduction chain, that 
is in Sp(10) model, the SU(5) Young labels agree with those given in the preceding section 
[[32]| . In the boson space, however, there is no Pauli-principle restriction. One immediately 
deduces that for iV < — , the correspondence between fermion states and boson states is 
one-to-one. For N > j , the most symmetric bosonic SU(5) representations are unphysical. 
In the fermionic SU(5) space, these representations are forbidden as the most antisymmetric 
conjugate representations are restricted by the dimensionality of the inert space. 

Turning to the SO(12) model, the SU(6) Young labels of the relevant U(15)dSU(6) 
reduction agree with the labels given in the preceding section for the fermionic active space 
In the boson space, there is again no Pauli-principle restriction. For N > ^, unphysical 



representations occur in the boson space. 

An alternative way to investigate the spurious states discussed above, is to calculate 
the overlap of generalized coherent states (|7.5f) . This overlap was calculated for the SO (8) 
and the Sp(6) models in Ref. Tailoring the discussion in Ref. [SO] to the coherent state 



overlap of the generalized Ginocchio model, one sees that it is possible to write 

(c, n\c, n) = ((c, n \c, Q )f /n ° (8.1) 

where ffo = i + | when K = and Qq = 2k + 1 when 1 = 0. Consequently, one expects 
spurious states to appear when N > where A" is the number of bosons. 

For the SO(12) model analyzed in terms of s-, d-, and g-bosons, the condition for the 
appearance of spurious states is thus A^ > ^ and for the Sp(10) model |TT| N > j. This 
should be compared (for the s- and d-boson models) with A^ > ^ for SO (8) and N > j for 
Sp(6). The condition for the SO(8) model in fact means that there are no spurious states 
when bosons are counted as in IBM, where hole bosons determine the boson number past 
mid-shell, implying that the number of bosons is always less or equal to a quarter of the 
shell size [p0| . Clearly this one-to-one correspondence between physical and ideal sci-boson 
states in SO(8) does not generalize to sdg-states in SO(12). 

In principle, one could use the U(15)dSU(5) and U(15)dSU(6) classification schemes of 
the ideal boson space in the Sp(10) and SO(12) models, respectively, and then restrict the 
physical space by throwing out the most symmetric unphysical SU(5) or SU(6) representa- 
tions. Algebraic techniques of the above mentioned classification chains are not, however, 
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developed in full details. More practical seem to be approaches in which d- and g-boson 
spaces are separated as e.g. in the U(15)dU(5)(S>U(9) chain. Then, however, physical and 
unphysical components mix in the boson basis states and only after diagonalization of the 
Hamiltonian they become separated. 

As for the actual identification of spurious states, these may be recognized by calcu- 
lating the matrix elements of a general SO(12) or Sp(10) operator between the eigenstates 
of whichever SO(12) or Sp(10) Hamiltonian is being used. This method is based on the 
structure |§ 

(^spurlOJVViys) = 0, (8.2a) 

(V ys |9 D |^ spur ) ^ 0. (8.2b) 

In other words, the Dyson boson image G D of any fermion operator O does not scatter 
outside the physical subspace when acting on a physical ket state. 

Another possibility to identify spurious states utilizes the application of the so-called 
7£-projection |30|, The 7£-projector acting on a boson state replaces each boson creation 
operator Bj M by the Dyson boson image -Rj M (|7TT|) . Because spurious bra-eigenstates are 
not contaminated by physical components (which is not true for the spurious ket-eigenstates) 
a spurious bra-eigenstate of the boson Hamiltonian is annihilated by the 7£-projection. A 
suitable way to detect spurious states in this way is to evaluate the norms of the 7£-projected 
bra-eigenstates of the Hamiltonian || |3"Ifl . The norms of spurious states are zero. To 
calculate the norms one needs the overlap of ^-projected basis states with non-projected 
basis states. This can be computed by iteration for different boson numbers from the relation 

K(N,L;I a ,I b ) = £ i— (AT/ 6 L||4||iV-lJ^') 

L'I'J't LL HJ 

(NI a L\\Rl\\N - U' a L')K{N — 1,1/; I' a , I' b ) (8.3) 

Here R\ is given by Eq. ( f7. l| a) and J should be chosen in such a way that the first re- 
duced matrix element in ( |3.3| ) is non-zero. The iteration starts from N = 2 and 7Z(N = 
l,L;I a ,Ib) = L5 a b- To decide whether an eigenstate \N, ip, L) is spurious or not we just 
calculate J2i b TZ>{N, L; I a , h)(N, (p, L\NI b L). This quantity is zero for each I a if the state is 
non-physical. 

In general, the method ( |3.2|) is simpler to implement. One has, however, to calculate 
quite a number of matrix elements to identify spurious states unambiguously. On the other 
hand, the ^-projection is somewhat more complicated to implement, but the identification 
is made from a single matrix element. From previous studies || |[J we should point out that 
the ^-projection is more reliable when one considers a slightly truncated boson system or 
one with slightly broken symmetry. 



Further application in the SO(12) model can be found in Ref. |JT| where it is demon- 
strated that a good fit to a spectrum from a phenomenological boson analysis does not 
necessarily guarantee the absence of large spurious components. 
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9 Quadrupole collective spectra in the SD and SDG 
models 



In this section we present examples of collective spectra which result from the SDG fermion 
pair algebraic models and compare them with the well-known SD solutions. We consider 
the N =6 major shell of normal-parity subshells, where one can construct the SD Sp(6) and 
SDG SO(12) models in the j=\, . . . ,f space (fi=21), and the SD SO(8) and SDG Sp(10) 
models in the J=|, ■ ■ ■ space (f2=20) (see Table |). 

All the calculations were done after the Dyson boson mapping had been performed in 
the ideal s, d or s, d, g space using modified and generalized interacting boson model codes. 
The ideal boson basis was classified according to U(5) chain for <i-bosons and U(9) chain 
for g-bosons. 

For the SD models we consider the Hamiltonian 

H = (1 - x) (G S + S + G 2 D + ■ D) + xk 2 Q ■ Q, (9.1) 

where S + S and D + -D are the monopole and quadrupole pairing interactions, respectively, 
and Q-Q is the quadrupole-quadrupole interaction. In notation of Sec. |3], the pairing and 
multipole operators are identified as S + =Fq Q , Z7^=F^q, and Qm=P2M with traditional 
normalization factors of g=£l and h=—\f2£L. Parameter x allows for a continuous change 
of the Hamiltonian between two limiting cases of a pure pairing-like interaction (x—0) and 
a pure quadrupole-quadrupole interaction (x=l). 

In Figs. and ^ we show the results for the SO(8) and Sp(6) models, respectively, 
obtained with the interaction strengths Gq=— 80 keV, G 2 =— 30 keV, and k 2 =— 25 keV. All 
calculations are for Np=8 fermions, i.e., for the boson number N=A smaller than ^ and § 
for the values f2=20 or 21 actually used. This means that in these cases all boson states 
are physical. Even if we use the same strengths parameters in both models, one should 
note that the operators in terms of which the SO(8) and Sp(6) Hamiltonians are defined 
are not the same. The phase spaces are also different, because the j—\ subshell is excluded 
in the SO(8) model for iVo=6. On the other hand, the size of the fermion space influences 
the results only through the parameter Q which appears in the boson mapping, and similar 
spectra can in fact be obtained for rather different values of Q. 

The left parts of Figs. and || show complete SO(8) and Sp(6) spectra for x=0 as 
functions of the spin /. In both cases one obtains characteristic vibrational patterns of 
quadrupole collective excitations. The ratios of excitation energies in the yrast bands are 
close to the vibrational limit of 2, 3, and 4 for 1=4, 6 and 8, respectively (see Table |TT|). 
For the SO (8) model, the Hamiltonian has exact SO (5) dynamical symmetry and many 
states of different spins are therefore degenerate. For Sp(6) this is not the case, but at x=0 
one can still clearly distinguish approximate multiplets corresponding to vibrational phonon 
excitations. 

The centre parts of Figs. [| and |3| show the dependence of spectra on the parameter x 
and exhibit the characteristic pattern of a phase transition to deformed structures. The 
energies of the first 2 + excitations decrease and the ratios of excitations in the yrast bands 
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increase. In the SO (8) case, the limiting situation of a pure quadrupole-quadrupole force 
(x=l) corresponds to the SO(6) dynamical symmetry and to the spectrum of a 7-unstable 
rotor with ratios of 2.5, 4.5, and 7 (Table fj). For Sp(6) the same limit gives the rotational 
structure resulting from the SU(3) dynamical symmetry. 

We may now compare the results of the SD algebraic models with those which use the 
SDG fermion pairs. For the N =Q major shell we have calculated the collective spectra in 
the SO(12) and Sp(10) models defined in their corresponding phase spaces listed in Table |. 
The boson number N=A is also here not larger than either ^ or j and therefore all boson 
states are physical. 

It turns out that the results for Hamiltonian ( |9.1| ) exhibit very low lying 4 + vibrational 
states which can be interpreted as hexadecapole excitations of the system. If one wants 
to discuss pure quadrupole structures one may push up the L=4 vibrations by adding an 
appropriately strong G + -G pairing interaction, where G^^F^ is the hexadecapole fermion 
pair creation operator. Therefore, in what follows we present results for the Hamiltonian 

H=(l-x) (G S + S + G 2 D + ■ D) + xk 2 Q ■ Q + G 4 G + ■ G (9.2) 

with coupling strengths Go, G 2 , and k 2 equal to previous values and G 4 being equal to 
180 keV and 250 keV for SO(12) and Sp(10), respectively. These values of G 4 ensure that 
the hexadecapole excitations do not appear below 3 MeV. 

In Figs. [| and ^] we show the results for the SO(12) and Sp(10) models using Hamilto- 
nian (|9.2|) . For x=0 we again obtain vibrational patterns with characteristic multiplets of 
quadrupole phonon excitations in which all degeneracies are lifted. In the SO(12) case, we 
see the well-isolated lowest hexadecapole state which at x=0 starts slightly above 4 MeV and 
at x=l goes down to about 3.2 MeV. When x increases from to 1 one obtains the phase 
transition to deformed spectra. The symplectic Sp(10) model now gives ratios of yrast ex- 
citations close to those of a 7-unstable rotor (Table ||), while the orthogonal SO (12) model 
ratios resemble those of a rotor. This situation is opposite to what one obtains in the SD 
models. 

10 Summary and conclusions 

In the present paper we have discussed algebraic collective models based on the concept of 
favored pairs of nucleons coupled to given angular momenta. First of all we have shown that 
the Ginocchio construction of such models has its roots in a simple assumption concerning 
the angular momentum symmetry. In fact, if one wants to replace (in a restricted single- 
particle space) the physical angular momentum by another J—l operator, and keep their 
matrix elements identical, this space necessarily splits into a set of r-p-shells of states for 
which two different angular momentum labels are simultaneously valid. This has the direct 
consequence of splitting the physical angular momentum into active and inert parts. 

Along these guidelines we have constructed generalized Ginocchio models for pairs of 
states with angular momenta varying between and an arbitrary value of J max - This leads 
to two families of algebraic models based on the SO(4p+2) and Sp(4p+2) symmetries for 
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half-integer and integer p, respectively, where J max =2[p]. We have also identified possible 
subgroup chains and representations as well as linked the model spaces to the nuclear shell 
model. 

Similarly as in other algebraic models, the structure of basic building blocks (here the 
structure of pairs) is fixed by the algebra closure conditions and takes precedence over 
dynamical considerations. We have analyzed this structure by discussing another set of 
building blocks, i.e., the monopole pair and the quadrupole single-particle operator. The 
latter has been compared with the mass quadrupole moment operator acting in one harmonic 
oscillator shell. This has demonstrated that the models based on the S and D fermion pairs 
lead to quadrupole operators rather different from the harmonic oscillator estimate. On 
the other hand, when the J=4, and then the J=6 fermion pairs are included the agreement 
improves. In this way the algebraic models may become closer to the traditional shell model 
description. 

By including higher and higher angular momentum pairs we have to consider larger and 
larger algebras, and the usual group theory techniques used to solve collective dynamical 
problems become more and more cumbersome. In the present paper we have demonstrated 
that the boson mapping technique can be a viable alternative. We have introduced Dyson 
boson mappings for the general class of algebraic fermion models considered here and then 
diagonalized Hamiltonians for two SDG models based on the SO(12) and Sp(10) symmetries. 
In fact the standard computer codes which have been constructed to solve the sdg-boson 
IBM models can be used with minor modifications. 

An obvious advantage of this procedure is the fact that one is not bound to the specific 
symmetry limits of the algebraic models and an arbitrary Hamiltonian expressed in terms 
of generators (not necessarily Casimir operators) can be rather easily diagonalized. 

As an example, we have applied this technique to Hamiltonians containing monopole 
and quadrupole pairing as well as quadrupole-quadrupole interactions. In particular, the 
SDG algebraic models based on the SO(12) and Sp(10) symmetries have been considered 
and compared with the well-known SD SO(8) and Sp(6) models. By changing the rela- 
tive interaction strengths between the pairing and quadrupole-quadrupole terms we have 
obtained phase transitions between anharmonic vibrations and rotations. Very good ro- 
tational structures have been obtained in the SO (12) model where there is no underlying 
SU(3) symmetry limit. We have found very low-lying hexadecapole vibrations present in 
both SDG models, which however could be pushed up to higher energies by using suitable 
hexadecapole pairing interaction. 

This work was supported by grants from the Foundation for Research Development of 
South Africa, the University of Stellenbosch, in part by the Polish State Committee for 
Scientific Research under Contract No. 20450 91 01, and by grant No. 202/93/2472 from 
the Grant Agency of the Czech Republic. 



19 



Appendix: Algebraic structure of the generalized Gi- 
nocchio model 



Fermion pairs of Eq. ( [3.3|) , which are the basic building blocks of the generalized Ginocchio 
model, can be presented in the following way: 



F JM = J o(2r + rt ^ 1 -il £ (pm p pm' p \JM) 



x (n + n + —( — 1) 2r n + n + \ (AD 

y rm r pm p rm r pmL \ ) rrn r pm p rm r pm p J i \ / 

where 

n + = (—l) r ~ mr n + (A 2) 

^ririrprrip V / r,— m r pm p " V / 

For different values of m r >0, the terms in Eq. ( |A.1| ) are independent of one another and 
the algebra comprised of the pair creation operators Fj~ M , their hermitian conjugates and 
all possible commutators is therefore a direct sum of algebras for all m r >0. 

For integer r we may further decouple the a and a fermions by using the unitary trans- 
formation: 

2(l + S mM ) (A3) 

( fl r,-ra r j)m p — ^■r,-m r prn p ) f° r m r — 0) 



^ rm r pm p 



/2(l + 5 mr0 ) 
the non-zen 

pair creation operator is a sum of 2r + 1 independent pairs for different m r values, i.e., 



where for m r =0 one takes the non-zero result. In terms of fermions a' „ m , the fermion 



JM 



oto i (pm p pm' p \JM)a' (A.4) 

V > J m r m p m' p 



Hence the algebra is the direct sum of 2r + 1 identical copies of the pair algebra for fermions 
occupying a single-p shell. Since this shell contains 2p + 1 states, the algebra in this case is 
so(2(2p+l)). 

For half-integer r the decoupling of the a and a fermions is not possible. However, in 
this case the terms in Eq. ( |A.1| ) are symmetric with respect to exchanging the indices m p 
and m' p and therefore they can be represented as pairs of bosons. In fact we may define 
these boson pairs as 



X J M = \ 0(ofj. ri ^ ^ (P m pP m 'p\ JM )^ mr p mp h lrn rPm U ( A - 5 ) 
V Z[ . Z ' l )m r >0 mpm p 



where &J mrpmp create bosons in the (2r + l)/2 copies of the single-p shell. Note that in the 
case of half-integer r there is no m r =0 term in Eq. ( |A.5|) . 



It is now easy to show that the algebra of the fermion pairs Fj~ M , Fjm and their commu- 
tators is identical to the algebra of the boson pairs Xj M , —Xjm, and their commutators. It 
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is important to note that one has to change the sign of the boson pair annihilation operator 
Xjm in order to obtain the identical algebras. It means that we use different real forms of 
both algebras and their representations are of course not equivalent. However, the angular 
momentum coupling does not depend on which real form is used. Hence the calculations 
can be facilitated by working in the angular-momentum-uncoupled representation and we 
may use the commutation relations of boson operators b^bl to calculate those of the fermion 
operators a+a+ + a+a+. For bosons occupying the 2p+ 1 states of the single-p shell the alge- 
bra of boson-pair operators is sp(2(2p+l)) and such is therefore the generalized Ginocchio 
algebra for half-integer r. 
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Table I: The SD, SDG, and SDGI fermion-pair algebraic models in normal-parity nuclear 
shells with radial quantum numbers iVo=3,. . . ,6, i.e., for particle numbers between 28 and 
184. The asterisks denote the cases studied in the FDSM fllQf . 
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Table II: Ratios of excitations Ri=(Ei+—E +)/(E2+—E +) in yrast bands of algebraic SD 
and SDG models in the vibrational (x—0) and rotational (x—1) limit. Results for the 
vibrational, rotational, and 7-unstable quadrupole model are also given for comparison. 



Model 


J = 4 1 = 6 1 = 8 I =10 


vibrational 
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3.16 6.47 11.19 22.40 
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2.5 4.5 7 10 



24 



Figure 1: Amplitudes defining the quadrupole operator Q, Eq. (|6.7|) , in the SD SO(8) and 
Sp(6), SDG SO(12) and Sp(10), and SDGI SO(16) and Sp(14) algebraic models constructed 
in the N =6 normal-parity major shell. The asterisks denote analogous amplitudes of the 
mass quadrupole moment in the Nq—6 harmonic-oscillator shell plotted in units of (r 2 ). 



Figure 2: Spectra of the SD SO(8) Hamiltonian for Np=8 fermions moving in the Q=20 
normal-parity major shell corresponding to the radial quantum number of iV =6. The 
interaction strengths of Hamiltonian ( |9.1| ) are Gq=— 80 keV, G2=— 30 keV, and K2=— 25 keV. 



Figure 3: Same as Fig. |2| for the SD Sp(6) model and Q=21. 



Figure 4: Spectra of the SDG SO(12) Hamiltonian for Np=8 fermions moving in the Q=21 
normal-parity major shell corresponding to the radial quantum number of Nq=6. The 
interaction strengths of Hamiltonian ( |9.2|) are Gq=— 80 keV, G2=— 30 keV, G4=180 keV, 
and K2=— 25 keV. 



Figure 5: Same as Fig. | for the SDG Sp(10) model with G 4 =250 keV and fi=20. 
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